Unconventional Neel and dimer orders in a spin-| frustrated ferromagnetic chain with 
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We study the ground-state phase diagram of a one-dimensional spin-^ easy-plane XXZ model 
with a ferromagnetic nearest-neighbor (NN) coupling Ji and a competing next-nearest-neighbor 
(NNN) antiferromagnetic coupling J2 in the parameter range < J2/IJ1 < 0.4. When J2/IJ1 < 
1/4, the model is in a Tomonaga-Luttinger liquid phase which is adiabatically connected to the 
critical phase of the XXZ model of J2 = 0. On the basis of the effective (sine-Gordon) theory and 
numerical analyses of low-lying energy levels of finite-size systems, we show that the NNN coupling 
induces phase transitions from the Tomonaga-Luttinger liquid to gapped phases with either Neel or 
dimer order. Interestingly, these two types of ordered phases appear alternately as the easy-plane 
anisotropy is changed towards the isotropic limit. The appearance of the antiferromagnetic (Neel) 
order in this model is remarkable, as it is strongly unfavored by both the easy-plane ferromagnetic 
NN coupling and antiferromagnetic NNN coupling in the classical-spin picture. We argue that 
emergent trimer degrees of freedom play a crucial role in the formation of the Neel order. 
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I. INTRODUCTION 

The search for novel orders arising from geometric 
frustration and quantum fluctuations in low-dimensional 
quantum spin systems has been a subject of intensive 
theoretical and experimental research. One-dimensional 
(ID) systems offer particularly ideal grounds for the- 
oretical studies, as powerful non-perturbativc methods 
are available (for a review, see, for example, Rcf. [l]). 
A prototypical example hosting a variety of fascinat- 
ing phenomena is the spin-^ frustrated chain with the 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
exchange interactions: 
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Here, Sj = (SJ , SJ , SJ) is a spin-i operator at the site 
j, and A is an XXZ anisotropy. The model has frustra- 
tion as long as the NNN coupling J 2 is antiferromagnetic, 
irrespective of the sign of J\ . 

Much of recent interest in the model (Q) has been fo- 
cused on the case when ferromagnetic (FM) J\ < and 
antiferromagnetic (AFM) J 2 > compete. The growing 
interest was triggered by experimental studies of various 
quasi-lD edge-sharing cupratcscHa which can be mod- 
elled by the Hamiltonian (|l|). Theoretical studies have 
uncovered a rich phase diagram of the model in a mag- 
netic field, which includes the vector chiral phase, spin 
nematic, and various other multipolar liquids .ErtL.3 Even 
without a magnetic field, the model (|l|) shows a vari- 
ety of phases depending on the values of the exchange 
anisotropy A and the frustration parameter JaJ \J\ \ . Ear- 
lier studies have discussed the diipep-prdcrOEj and the 
stability of the vector chiral ordertlM in the easy-plane 



case < A < 1, and the partial ferromagnetismE! in the 
easy-axis case A > 1. In this paper we show that the 
phase diagram of the J\-3i chain with easy-plane XXZ 
anisotropy, Eq. ([j]), exhibits unexpectedly complicated 
and interesting phase structure, where Neel and dimer 
ordered phases alternate more than twice, near a FM 
critical point (J2/\Ji\, A) = (0.25,1). Our analysis is 
based on the combination of bosonization and numerical 
analysis of finite-size energy spectrum, which has been 
proven successful in the study of many ID quantum spin 
models including the model (Q) with the antiferromag- 
netic NN coupling. 

The ground-state phase diagram of the anisotropic J\- 
J 2 spin chain (|j]) is well understood in the antiferro- 
magnetic case, J\ > and J 2 > (see, e.g., Refs. p^ - 
prj| ). For small J%j J\ and with easy-plane anisotropy 
A < 1, the model is in a gapiess phase which is de- 
scribed as a Tomonaga-Luttinger liquid (TLL). The gap- 
less TLL phase has instabilities towards antiferromag- 
netic (Neel) order and dimcrizationpdis exemplified by 
the exact singlet dimer ground statdlj for J%j J\ = 1/2 
and A > —1/2). Both Neel and dimer ordered phases 
are gapped. The Neel order is induced by easy-axis 
anisotropy A > 1, while the dimer order appears when 
J2/J1 i s MfSf r than some critical value (rj 0.24 at A = 1). 
HaldaneEJEll showed that the quantum phase transitions 
to these ordered phases can be understood in a unified 
way within the sine-Gordon (SG) model which is ob- 
tained by bosonizing the Hamiltonian (jl|). The two types 
of orders arise as a bosonic field is locked at two distinct 
values by the cosine potential coming from the Umklapp 
scattering of Jordan- Wigner fermions. It is also worth 
noting that, for large J 2 / J\, a vector chiral ordered phase 
with gapiess excitations appears for a certain range of 
easy-plane anisotropyj^jl 2 ^ while thcrdimcr ordered phase 
persists in the isotropic case A = l.cil 
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In this-j-jpaper, we apply the SG formalism of 
HaldaneEltj to the frustrated ferromagnetic spin chain 
(Q). Our starting point is the FM Ji-only chain with 
easy-plane anisotropy < A < 1. In this case the system 
is a TLL whose properties are understood in great detail 
from the Bethe ansatz and the bosonization. In particu- 
lar, thejaai3,meters in the effective SG theory are known 
exactly.E3c3 According to these exact results, rather in- 
terestingly, there occurs a scries of sign changes in the 
coupling of the cosine potential, as a function of the 
anisotropy A. In the J\ XXZ chain, the cosine poten- 
tial is irrelevant in the renormalization-group (RG) sense 
and does not make any significant impact on the long- 
distance properties. However, it is made relevant by large 
enough 3 2 coupling, and as a result we observe alternat- 
ing appearance of Neel and dimer ordered phases, as the 
exchange anisotropy is changed towards the FM critical 
point (Ja/IJiJ, A) = (0.25,1). We determine the phase 
boundaries accurately using finite-size scaling analysis 
of discrete ^energy levels, the so-called level spectroscopy 
method, 000 which combines the SG theory with the 
numerical exact diagonalization. The appearance of the 
antifcrromagnetic spin ordering in the z direction is par- 
ticularly counterintuitive; on the classical level, such spin 
configuration is strongly unfavored by the easy-plane FM 
Ji coupling as well as by the AF J 2 coupling. We ana- 
lyze the correlations and the reduced density matrices 
of the ground state in the Neel phase using the infinite 
time-evolving block decimation (iTEBD) algorithm,E3 a 
numerical method which can directly address physical 
quantities in the thermodynamic limit. We will then ar- 
gue that emergent trimer degrees of freedom play a cru- 
cial role in the formation of the Neel ordered state. 

The present study focuses on the region with J2/I Ji| < 
0.4. The analysis of the region with larger J 2 /\Jii-i& also 
done by iTEBD and will be reported elsewhere£2l It is 
found that for J2/\Ji\ > 1/4, the vector chiral ordered 
phase with gapless excitations is robust up to a very weak 
anisotropy Awl and appears in a wide range of the 
parameter space. This is in contrast to the antifcrromag- 
netic J1-J2 model, J«here the vector chiral ordered phase 
has been idcntificcO only for large J2/Ji(> 1-2). 



II. FORMALISM 

A. Effective field theory 

We first consider the case J2 = 0, i.e., the easy-plane 
XXZ chain with J\ < 0. We parametrize the exchange 
anisotropy as 



A = cos(7r?7), < 77 < i. 



(2) 



In this parameter range, the model is in the gapless TLL 
phase. In the bosonization formalism (for a review, see, 
e.g., Rcfs. |32],p3|), its low-energy effective theory is the 



quantum SG theory defined by the Hamiltonian density 
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where the bosonic field (j> and its dual counterpart 9 
obey the commutation relation [4>(x),9(y)] = — (i/2)[l + 
sgn(a; — y)] , and the lattice constant is set to unity. The 
TLL parameter K and the spin velocity v are given by 
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The cosine potential in Eq. (J3|) has scaling dimension 4K 
and is irrelevant in the XXZ chain with the exchange 
anisotropy (||). We keep the cosine term because it will 
drive the phase transitions from the TLL to the Neel and 
dimer ordered phases in the presence of the frustrated 
NNN coupling J 2 > 0. When J 2 = 0, tl^xact value of 
the coupling constant A is known to becJ'Ej'tfJ 
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where the normalization condition on the zero- 
temperature correlator of vertex operators, 

(e^'e-^l*'') = \x- x'\- K fl^, \x -x'\ > 1, (6) 

is assumed. Notice that Ao vanishes at 77 = 1/n, i.e., 

A = cos (71-/71) , 71 = 3,4,.... (7) 

At these points, the XXZ model is invariant under the 
action of the loop algebra sl 2 , leadins-tp-some non-trivial 
degeneracies in the energy spectrumO'CJ We will observe 
one example of such degeneracies in our numerical result 
later. 

The spin operators of the original lattice model (|l|) are 
expressed ast3 



1 d(j> 
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where a, 60, and 61 are non- universal constantsB^S Sim- 
ilarly, the staggered part of the NN exchange energy is 
given by 



Sj ■ Sj-! - Sj ■ S j+ i = c(-l) j sm(V4^<f>) 



(10) 



where c is another non-universal constant. 

Including the NNN coupling J 2 > changes the pa- 
rameters (v, A, and K) in the effective theory (|3). A per- 
turbative calculation using (||) and ([)]) yields the renor- 
malized TLL parameter for J2/IJ1I 1, 
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FIG. 1: RG flow diagram of the sine-Gordon theory (|3(). 
The thickr-solid lines indicate the Kosterliz-Thouless-type 
transitional^ between the TLL and the two ordered phases. 
The dotted line indicates the Gaussian transition between the 
two ordered phases. 

where Kq and vq are the values of K and v at Ji — given 
in Eq. (||). This indicates that K decreases as J2/\Ji\ 
increases. As long as K is larger than 1/2, the cosine 
term is irrelevant and the system is in the TLL phase. 
As J2/\Ji \ is further increased, the cosine term in Eq. (||) 
will eventually become relevant (AK < 2). The positive 
and negative values of A will then trigger the Neel and 
dimer orders, respectively, as one can see from Eqs. (j^) 
and ([lO]). Wc analyze such ordering transitions through 
a combination of the RG analysis of the SG model (|3|) 
and the finite-size spectra obtained from numerical exact 
diagonalization. 

B. Renormalization group flow 

Wc discuss the RG flows of the SG theory (|) by vary- 
ing a short-distance cutoff a. To this end, it is convenient 
to parametrize 

a = a e, A =- + —-, A= — — , (12) 
2 4 or 

where ao is the initial cutoff. Under an infinitesimal 
change I — > I + dl, the dimensionless parameters y(l) and 
y ^(l) change according to the RG equations 

dy{l) 2 /n C M0 m m mi 

-^- = -2/0(0: —^—=-^(02/(0- ( 13 ) 

The RG flow diagram is drawn in Fig. |l|. When < 
y < | 1, the cosine potential finally vanishes, and the 
system is described as the TLL. The transitions from the 
TLL to the Neel and dimer phases occur along the half 
lines y$ = +y and y^ = —y with y > (thick solid 



FIG. 2: (color online) Energy difference (L/v)(Eo — -En) 
in the unfrustrated case J2 = 0. The arrows indicate the 
"special" points in Eq. (Q) with n = 3, 4, . . . , where A changes 
the sign. 

lines), respectively. The transition between the two or- 
dered phases occurs along the Gaussian fixed line y^ = 
with y < (dotted line). 



C. Level spectroscopy 

On the basis of the RG flow diagij^an of the SG model in 
Fig. [I], Okamoto and NomuraE3n3c3 developed a simple 
numerical method which allows to determine the phase 
transition lines precisely The key idea of this method 
is to look at the lowest excited states in the three dif- 
ferent phases. Under the periodic boundary condition, 
the eigenstates of a finite-size system of length L arc 
labeled by the magnetization S z = J2j &j > the wave 
number q = 2irk/L (k £ Z), the (bond-centered) par- 
ity P = ±1, and the spin reversal T = ±1. In the pa- 
rameter region of our interest, the ground state (with 
energy E ) of the model (jl]) is in the sector (S z = 0, q = 
0, P = +1, T = +1) when L is even. In the TLL phase, 
the first excited state (with energy E$) is in the sector 
(S z = ±1, q = 0, P = +1). In the Neel and dimer phases, 
pscudo ground states (with energies and Ed) appear 
in the sectors (S z = 0,q = tt,P = — 1,T = —1) and 
(S z = 0, q = 7r, P = +1,T = +1), respectively, reflecting 
the two-fold ground-state degeneracy in the thermody- 
namic limit . One can expect that the transitions between 
the three phases can be detected by observing crossing 
of these energy levels. 

This expectation can be justified by a detailed aHalysis 
of the finite-size spectra of the sine-Gordon theoryBj See 
Appendix for details. Here we briefly mention the basic 
idea. The excitation energy AE§ := Eg — Eo is associated 
with the operator ~ e ±l V^ e ( x )^ The excitation ener- 
gies Ai?N := -En — Eq and A-Ed := -Ed — Eq are associ- 
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ated with the order operators (— l) 3 S* ~ cos(V47raV) and 



(—l) J (Sj~Sj +1 + h.c.) ~ sin(-\/47r</>), respectively. From 
the correspondence between the spectra and the opera- 
tors, one obtains 



AE n = 



2™ fi_m\ 

L \2 4 J ' 
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(14a) 
(14b) 
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where Z is related to the system size L by e' = L/2ir. 
One can see that the crossing of Eg and En occurs when 
y(l) = y<f,(l), corresponding to the TL-Neel transition in 
Fig. [lj. Similarly, Eg = Ed corresponds to the TL-dimer 
transition. Using the formula L(Ed — £n) = 2itvy<f > (l) 1 
one can observe the running coupling constant, in partic- 
ular, its sign change. In Fig. ^, (L / v) (Ed — -En) is plotted 
for the J i chain (J2 = 0). For a system of L spins, we 
observe zeros of y<t>(l) at "special" points in Eq. ([?]) up 
to n = L/2. In the faiit L — > 00, the energy difference 
asymptotically obeysEJ 



'■(Ed-Es) =27tA 



2-2/t, 



(15) 



This formula is also plotted for L = 24 (red solid line in 
Fig. ||), which agrees well with the numerical result for 
A < 0.8. For larger A, better agreement is expected to 
be seen at larger L. 



III. PHASE DIAGRAM 

The phase diagram is determined by the level spec- 
troscopy and the iTEBD analysis as presented in Fig. ||. 
The blue "+" symbols show Gaussian fixed points deter- 
mined from the condition E^ = Ejj . The black open cir- 
cles interpolated by lines show the transition points from 
the TLL phase to the Neel and dimer phases, which are 
determined by Es = £"n and Eg = Ed , respectively. This 
phase boundary between the TLL phase and the gapped 
ordered phases is smoothly connected to the ferromag- 
netic transition poinica J2/IJ1I = 1/4 in the isotropic 
case A = 1, at which the ground states are known to be 
highly degenerated The Gaussian fixed lines with A = 
start from the "special" points given in Eq. (^) at J 2 = 0, 
and are therefore labeled by n = 3, 4, . . . for convenience. 
Some of them continue even after the cosine potential 
becomes relevant for J2/I Ji| > 1/4. From calculation of 
L = 24 system, we have identified five lines (correspond- 
ing to n = 3, . . . , 7) extending into the parameter region 
h/\J\\ > 1/4 and A < 0.995. Across these lines, there 
occur successive transitions between the Neel and dimer 
phases as the anisotropy parameter A is changed toward 
unity. However, for 0.995 < A < 1, Lanczos diagonal- 
ization does not converge well, due to highly degenerate 
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FIG. 3: (color online) Phase diagram of the model (fil). The 
level spectroscopy has been performed for L = 24. We have 
observed that finite-size effects are quite small and can be ne- 
glected. The inset shows a zoom around the highly-degenerate 
point (J 2 /|Ji|,A) = (1/4,1). 



nature around (J2/|Ji|,A) = (1/4,1). For this reason 
we have not been able to find the transition lines beyond 
n = 7. 

Inside the Neel and dimer phases, there is a Lifshitz 
line (pink "x" symbols), where the short-range spin cor- 
relation in the xy plane changes its character from com- 
mensurate to incommensurate (indicated by C and IC 
111 Fig. |). This line was determined by observing the 
peak position of tho-equal-time spin structure factor cal- 
culated by iTEBDfiil The level spectroscopy presented 
above is valid on the left side of this line. By further 
increasing J2/IJ1I, the system enters the vector chiral 
ordered phase in which the vector product of neighbor- 



ing spins, Kj d+1 



((Sj x Sj + i) z ), is long-range or- 



dered. The boundaries above and below the vector chi- 
ral phase (black diamond symbols) were determined with 
reasonable accuracy by observing the rapid increase of 
the order parameter Kj j+i calculated by iTEBD.tJ For 
0.9 < A < 1 and 0.25' < J 2 /|Ji| < 0.35, we have not 
been able to locate the transition lines accurately be- 
cause of a relatively poor convergence of iTEBD around 
the highly degenerate point (J2/|Ji|,A) = (1/4,1). How- 
ever, we expect that both the lines should continue to the 
point (J2/|Ji|,A) = (1/4,1). The vector chiral phase 
is smoothly connected to that found in the. Icuatrated 
ferromagnetic chain in a low magnetic fieldOEj'Eil Simi- 
lar vector chiral phases have also been found previously 
in the ant iferij omagnetic J1-J2 modgL^with easy-plane 
anisotropyOcj or in a magnetic field.l3cJL3 

The nature of the dimer phase at A < 0.7 is easy to 
understand, as this phase extends to the XY point A = 
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FIG. 4: (color online) The dimer order parameters and 
— 2Df 2 3 (left axis) and the local magnetization (Sf) (right 
axis), along the line connecting (J2/|Ji|,A) = (0.336,0.65) 
and (0.264, 0.95), which approximately traces the Lifshitz line. 
The order parameters take on relatively large values along 
this line. The calculation was done by the iTEBD with the 
Schmidt rank x = 300. The vertical lines indicate the phase 
boundaries determined by the level spectroscopy. 



0, where the sign of 3\ can be reversed by performing 
the 7T rotations around the z axis of the spins on every 
second sites. From the fact that the ground state at the 
Majumdar-Ghosh point J 2 = Ji/%-> and A = (and 
more generally for all A > — 1/2)E3 is given exactly by 
the product of singlet dimers (|t4-) — I4t))/V2, one finds, 
through the above 7r-rotation transformation, that the 
ground states at J 2 = — J%/2 > and A = are given by 
the dimer states whose dimer unit is now replaced by the 
triplet state (| tl) + I 4-"t))/v / 2-H We expect that such a 
"triplet" dimer nature should survive away from the XY 
point and define the order parameter of the dimer phase. 
The nature of the Neel phase in 0.7 < A < 0.9 is more 
elusive and will be discussed in detail in the next section. 

So far our argument for the emergence of the Neel and 
dimer phases has been based on the effective SG the- 
ory. To confirm their existence in an independent and 
unbiased way, we have calculated the ordcii-paramctcrs 
of these phases using the iTEBD algorithmic which can 
address physical quantities in the thermodynamic limit 
directly through the use of the matrix product represen- 
tation of (ground) states. When this algorithm is per- 
formed in an ordered phase, a variational state finally 
converges to a symmetry-broken state with a finite spon- 
taneous order parameter (if it is allowed by the periodic- 
ity of the matrix product state). Hence, the Neel phase 
is detected by a local magnetization (S*), and the dimer 
phase by dimer order parameters: 
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These order parameters show alternating signs along the 
spin chain. We choose the site labellings in such a way 
that (Sf) > and D\ 23 < 0. Wc calculated these order 
parameters for the parameter points (J2/I Ji\i A) varying 
approximately along the Lifshitz line; see Fig. EL In this 
figure we observe two transitions between the dimer and 
Neel phases, in agreement with the level spectroscopy 
analysis. In the dimer phases, D^ 23 and Df 23 show op- 
posite signs in accord with the "triplet" nature of the 



dimers. [Notice that D\ 



■>■<> 

23 



-2Df 23 = 1/2 in the exact 



"triplet" dimer state at (J 2 /|Ji|, A) = (1/2,0).] 

The sudden changes of the order parameters in Fig. Ej 
seem to indicate first-order nature of the transitions be- 
tween Neel and dimer phases for this parameter regime. 
On the contrary, a continuous transition of Gaussian type 
is expected in the SG model (Fig. 0). This discrepancy 
may be reconciled by considering the effect of a higher- 
frequency cosine potential cos(S^/n(f>), which was ignored 
in the SG theory but is allowed by symmetry, and can 
become relevant deep inside the ordered phases. With 
a negative coefficient, the potential cos(8v / 7T(/)) has four 
minima corresponding to the Neel and dimer ordcrings 
with doubly-degenerate ground states each. Different 
signs of A select different types of orderings, and thus the 
first-order transition at A = separates the two phases. 

We finally note that, in the small region near the up- 
per right corner of Fig. a weak dimcrrzation with the 
same sign in D* 23 and Df 23 was observed. Ell In particular, 
D*2 3 = 2D\ 23 in the isotropic case A = 1. This result 
indicates that this region is characterized as a dimer or- 
dered phase having a distinct nature from the "triplet" 
dimer phase in A < 0.7. 
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IV. NATURE OF THE NEEL PHASE 

The appearance of the Neel phase is a natural con- 
sequence of the SG theory (||) with A > 0, but is quite 
counterintuitive in the presence of the easy-plane nearest- 
neighbor ferromagnetic coupling. In Fig. Ej, the Neel or- 
der parameter is one order of magnitude smaller than 
the value (Sf) = 1/2 of the pure Neel state 1 . . . ), 
which indicates that the ground state is not well approx- 
imated by such a simple product state. Here we ana- 
lyze the nature of this phase in detail using iTEBD with 
Schmidt rank \ = 200. We first analyze the spin cor- 
relations; see Fig. ||. We observe that the Neel charac- 
ter of the correlations appears only at a relatively long- 
distance scale r > 6. At a short-distance scale, (SfSf +r ) 
takes larger amplitudes than (S*S* +r ), reflecting the 
easy-plane anisotropy. Furthermore, we note that (SfS 3 ) 
takes relatively large negative value, which is unusual in 
a Neel phase. 

Such unusual correlations in a short-distance scale can 
be understood in terms of trimer degrees of freedom. 
First, we consider a three-spin model shown in Fig. [3](a). 
In the isotropic case A = 1, the ground state is a quadru- 
(17) plot |q M ) (with S? ot = fi = ±§,±f) when J 2 /\Ji\ < 1/2. 
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FIG. 5: (color online) Spin correlation functions (SfS^+i) 
(left axis) and (S'fS'jr+i) (right axis) for J2/IJ1I = 0.31 and 
A = 0.74, where the Neel order parameter is relatively large: 
(ST) = -(Si) = 0.04432. 
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FIG. 6: (a): Three-spin model, (b), (c), (d): Neel order from 
the trimer picture. The arrows show the z-component of the 
local magnetization. 



Because of the reflection symmetry around the site 2, 
the excited states are classified into symmetric and an- 
tisymmetric doublets, and {|g?^)} (with /i = ±i). 

" \d. 



When the anisotropy A is introduced, |<Z±i) and 
are mixed to form new eigenstates |<?±i) and \d±i 
eigenstates are summarized as follows: 
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with 7 = i(VA2 + 8- A) and A = A+(J 2 /| J X |)(A-1). 
Other eigenstates with /1 = —1/2 or —3/2 arc ob- 
tained by applying the spin reversal to the above. The 
cigenencrgies are plotted in Fig. |t]. Under the easy- 
plane anisotropy A < 1, the ground states are given by 
\q±k (7)) (with 7 > 1), in which 5J's correlate ferromag- 
netically and 's antiferromagnetically for any spin pair 

(if 7 < \/2). This property is in accord with the anoma- 
lous spin correlations in Fig. || for r = 1,2. Note, how- 
ever, that the local magnetization (Sj) is almost uniform 
in \q±i (7)) as shown in Fig. ^](a), and we do not observe 
Neel ordering pattern at this level. 

Now we construct a Neel ordered state in a ID chain 
as follows. (The procedure is analogous to the construc- 
tion of the valence bond solic 
Kcnnedy-Licb-Tasaki model, 



sround state in the Afflcck- 



and more generalbi-jto the 
construction of a projected entangled-pair state.cS) The 



ID model ([j]) can be viewed as the model on the zigzag 
ladder, where J\ and J2 give the inter- and intra-chain 
couplings respectively. We place a trimer state (7)) 
(7)) resp.) on every up (down resp.) triangle in the 
zigzag ladder, as shown in Fig. ^](b). At this point, ev- 
ery site is shared by three neighboring trimers. To define 
the state in the original Hilbert space, we bind the three 
spins at every site together to form a single spin- 1/2, as 
indicated by an ellipse in Fig. |^(c). This is done by a pro- 
jection operator 1 7) (ip + | + 1 4.} \ , where | ip± ) are defined 
for the three spins in an ellipse. Wc require that (i) \t/j±) 
have the reflection symmetry about the central spin and 
that (ii) they have sufficient overlap with |f|t) an d 14-t-J-) 1 
which are expected to have large weights in the state be- 
fore projection as can be seen in Fig. g(c). The states sat- 
isfying (i) can be written in the same way as in Eq. (|l8|a)- 
(p"8|c). Considering also (ii), we set \i()±) = |g? ± i(7')) with 
7' appropriately tuned. The obtained state after the pro- 
jection has a Neel ordering pattern as shown in Fig. ||(d). 
We expect that this should give an approximate ansatz 
(with two variational parameters, 7 and 7') of the Neel 
ground state found in the present model. 

To support the trimer picture, we calculated the re- 
duced density matrix ^123 for neighboring three spins 
(1,2,3) from the ground state obtained by the iTEBD. 
It was done for the same parameter point as in Fig. |^. 
Because of the reflection symmetry about the site 2, /5123 
is diagonalizcd in the form 

P123 H<?+§)a+(9+§l + |5+i(7+))M£fi(7+)l 

+ |d + i( 7+ )}c + (d + i(7+)| + \d' H )e + (d' H \ (19) 
+ (" + " -> "-"), 
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FIG. 7: Spectrum of the three-spin model in Fie, H(a) for 
■/2/IJ1I = 0.31. The eigenstates are given in Eqs. (118^7 



where 7+ = 1.225 and 7_ = 1.250 arc not far from 7 = 
1.122 obtained in the three-spin problem. Large weights 
are found in b+ = 0.5803, &_ = 0.2633, e_ = 0.0999 
and a_ = 0.0472, and the others are orders of magnitude 
smaller. We see that the trimer states |<?±i(7±)) occupy 
about 84% of the total weight. 

We also applied the systematic method for determining 
the order parameter proposed in Ref. [i9[ In this method, 
we use the pair of symmetry-broken ground states ^1,2) 
- one is obtained directly by iTEBD and the other by 
applying the spin reversal to it. For a subregion A of the 
system, we construct the reduced density matrices p^ 2 
from these states, and measure the mutual distance: 



(20) 



where Xj's are the eigenvalues of p\ — p\. This measure 
takes a value from zero to two and shows to what extent 
the two states are distinguished in the concerned region. 
[See Rcfs. |49}|50| for details.] For one- and two-site regions, 
we obtain the same value diff = 4(5*) = 0.1773. The 
three-site region A = {1, 2, 3} gives a much larger value 
diff = 0.9128. These results indicate that, although the 
symmetry breaking can be detected through a single-spin 
operator as expected for a Neel phase, it is much more 
visible in terms of a three-spin operator. This also sup- 
ports the trimer picture for the Neel state. 

Finally, we note that a Neel ordering in the presence of 
easy-plane anisotropy has been found in .aa antifcrromag- 
netic model with explicit trimcrizationH This example 
also suggests that the formation of trimers can change 
the nature of the anisotropy of the system. 



CONCLUSIONS 



We have analyzed the phase diagram of the spin- 
\ frustrated ferromagnetic chain ([j]), starting from the 
bosonization picture in the unfrustrated case J2 = 0. A 
crucial observation is that the effective sine- Gordon the- 
ory shows a cascade of sign changes in the coupling of the 
cosine potential when changing the anisotropy A. When 
the cosine potential is made relevant by the competing J 2 
exchange interaction, it gives rise to alternate appearance 
of Neel and dimer orders. Although the Neel phase has 
a long-ranged staggered spin ordering pattern, we have 
argued that this ordering should rather be interpreted as 
the ordering of the trimer degrees of freedom. 

The successive sign changes in the coupling of the co- 
sine term in the effective theory of the FM J\ XXZ chain 
is explicitly seen in the exact expression (|^) obtained in 
Refs. ^7],|8| Our study has shown for the first time that 
it has a dramatic physical consequence, the alternation 
of the Neel and dimer ordered phases. It will be inter- 
esting to further examine the effect of the sign change in 
the presence of other types of perturbations. 
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Appendix: Finite-size spectra in the sine-Gordon 
theory 

In this appendix, following CardyS and Bortz et a/.,@ 
we present a simple derivation of the excitation energies 
( |l4| ) based on the perturbation theory from the Gaus- 
sian model. This is complementary to the derivation by 
NomuraEi using the correlation functions. 

In a finite-size system of length L, the effective Hamil- 
tonian is given by 



H =H< 



Gauss 







C dx < 




C 2 . 



dey 1 

\dx ) K \dx 



2tt 



where the lattice spacing is set to unity. 

Let us start from the Gaussian part -ffcauss 
to rescale the fields as <f> = 4>/VK and 9 
simplify the Hamiltonian: 



(Al) 



It is useful 

= Jko to 



H, 



Gauss — 



dx 



dO 

dx 



dx 



(A.2) 



8 



Assuming the periodic boundary condition, we perform 
the mode expansions of the bosonic fields: 



3 -M/2A 



$(x) = $ + Qy + V ^== (e ik " x a n + e~ ik ^n\ ' 



(A.3) 



e(a:) -Go + Q^ 



e -|n|/2A 

(A.4) 

with k n = 2im/L. Here, a n with n > (n < 0) rep- 
resents a right (left) moving mode. We have introduced 
an ultraviolet cutoff A, which will be determined later. 
The operators in the expansions obey the commutation 
relations 



[*o,e ] = -*, [Q,e ] = i, [Q,*o] = *, 



a n ,a 



(A.5) 



Using the mode expansions, the Hamiltonian is diagonal- 
ized as 



_ 2nv 

n Gauss — 



^(Q 2 + Q 2 ) + EH a « a » 



(A.6) 



where A is taken to infinity. The ground state |0) is 
defined by a„|0) = Q|0) = Q|0) =0 and its energy is set 
to zero. The compactification radii R and R for $ and 
O are given by 



2K 



1 



2-kR 1 



(A.7) 



2ttR 2 

Then Q and Q are quantized as 

Q = S z /R, Q = m/R, with S z ,meZ. (A.8) 

Here, the integer S z coincides with the total magnetiza- 
tion. The entire spectrum is then given by 



AE (S Z , m,{m n }) 



2ttv 



1 f(S z ) 2 m 2 



4tt V R 2 R 2 



lit \ sr-^ 

m + E 



n\m n 



(A.9) 



with m n > integer, and the corresponding eigenstates 
are 

\S z ,m,{m n }) = e iS°e /R+im* /R JJ( a t) m »|o). 

(A.10) 

The excitation AE$ is associated with e ±l6o ^ R and is 
thus given by 

A£ S = A^ ± 1,0,0) = ^ * =^ * (All) 

L AttR 2 L 2 + y 



with y = AK — 2. This coincides with the expression 
(14a) when \y\ <C 1. On the other hand, Ai? N / D are 
associated with cos($o/i?) and sin($o/-R) respectively. 
In the Gaussian model, these energies are degenerate and 
are given by 



A£ N/D = A£(0,±1,0) 



2nv 1 
L 4ttR 2 



2ttv 



K. (A.12) 



So far, we have ignored the cutoff A. It is determined 
so that the two-point correlation function of vertex op- 
erators satisfies the simple normalization condition 



( | e v*(*) e -v*(*')| ) = 



1 



(A.13) 



for \ x — x' \ <C L. Using the expansion (A.3), the l.h.s. of 
Eq. ( A.13 ) is evaluated as 



cxp [-n 2 (0\<f>{x)<$>(x')\0}] 



exp 



-ft > : | | 1 - e 



27ri\n\(x — x f ) I L 



1 . 2tt(x-x') 

1 — e A 1 



1 - 



2A; 



-Ai 2 /2ir 



V/27T 



(A.14) 



Hence, we set A = L/27T to satisfy Eq. ( A.13| ). 

Now we introduce i? C os as a perturbation. It splits 
the degeneracy between AE(Q, ±1,0). We calculate the 
matrix element: 



C := (0,-l,0|# CO s|0,l,0) 



'2tt 



dx (Ole^/^cos 



,/R 



|0). 

(A.15) 



Using the expansion (A3), the integrand is evaluated as 



2$(x) 



|0) 



| g i$o/ii e -2i*( ; i ! )/fl e i$o/ii|o) 



exp 



2 
1 

2< 

r xp (-^ 



2i 
R 



E 

oo 

E 



-H/2A 



,2^nx/L , jj^ 



|0) 



-A 



-4if 



(A.16) 



In first-order perturbation theory the new eigenstates are 
given by the linear combinations 



1 



V2 



(|0,1,0)±|0,-1,0)), 



(A.17) 
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and the corresponding cigcnencrgies are 
2ttv 



N/D 



-K±C 



2irv 



(A.18) 



Equations (A. 11) and (A.18) are non-perturbative in y 
and perturbative in A. 

So far, we have used the bare coupling constants, y and 
A. Now we relate the above results to the running cou- 
pling constants, y(l) and y<j,{l), focusing on the vicinity of 
the multicritical point (y,y<f,) = (0,0). Both the cosine 
potential and the change in the Luttinger parameter y 
are marginal perturbations at this point. In general, the 
excitation energy AE n associated with an operator O n 
is related to the inverse of the correlation length £„ of 
the operator. Under a global scale transformation by a 
factor e , it scales as 

AE n cx C X (2/, V* L- 1 ) = e- l ^- l {y{l),y^l), ZrV). 

(A.19) 



Setting e l = L/2w, this is given by a universal function 
in terms of y(l) and y^it) : 



AE n = ^ n (y(l) >y<p (l)). (A.20) 



Namely, LAE n depends on L only through y(l) and y$(l). 

We compare this with Eqs. ( A. 11 ) and ( A.18 ). As 
shown in Fig. [l], for y > \y^>\, y(l) finally converges to 
a constant. If the bare value of y is sufficiently close to 
the final constant, the RG equations are solved as 



2/(0 = y, y</>(i) = V4> e 



-yl 



(A.21) 



Thus , und er the correspondence e l = L/2n, Eqs. ( A. 11 ) 
and ( A. 18] ) are expressed by Eqs. which s how a lin- 
ear dependence on y(l) and y$(l). Equation ( A.20|) im- 
plies that we can extend the range of validity of Eqs. ( p4| ) 
to all |y(J)|,M0l«l- 
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